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Abstract. The Lorentz gas is a model for a cloud of point particles (electrons) 
in a distribution of scatterers in space. The scatterers are often assumed to be 
spherical with a fixed diameter d, and the point particles move with constant 
velocity between the scatterers, and are specularly reflected when hitting a scat- 
tered There is no interaction between point particles. An interesting question 
concerns the distribution of free path lengths, i.e. the distance a point particle 
moves between the scattering events, and how this distribution scales with scat- 
terer diameter, scatterer density and the distribution of the scatterers. It is by now 
well known that in the so-called Boltzmann-Grad limit, a Poisson distribution of 
scatters leads to an exponential distribution of free path lengths, whereas if the 
scatterer distribution is periodic, the distribution of free path behaves asymptoti- 
cally like a Cauchy distribution. 

This paper considers the case when the scatters are distributed on a quasi 
crystal, i.e. non periodically, but with a long range order. Simulations of a 
one-dimensional model are presented, showing that the quasi crystal behaves 
very much like a periodic crystal, and in particular, the distribution of free path 
lengths is not exponential. 



1. Introduction 

The Lorentz gas is a mathematical model for the motion of (point) particles in e.g. a 
crystal, consisting of spherical, elastic scatterers of radius a with centers at a fixed set of 
points T e C W 1 . A point particle moves in straight lines between the obstacles, on which 
it is specularly reflected. At least two very different scatterer distributions, T £ have been 
studied thoroughly: the standard lattice C C W 1 , with interstitial distance e( n_1 )/ n , or a 
random distribution, where T £ is Poisson distributed with intensity e - '" -1 ). In this paper 
we are mainly concerned with the so called Boltzmann-Grad limit of this system, and that 
corresponds to setting the scatter radius a = e and letting e 0. Very loosely speaking, in 
this scaling the mean free path of the point particles remain of order one as e — > 0. Consider 
now a density f E (x,v,t] of point particles moving between the scatterers. Gallavotti [9| 
(see also ifTUl ) considered the random obstacle distribution, and proved that when e — > 
0, f e (x,v,t) converges to a function that solves a linear Boltzmann equation. On the 
other hand, Golse ifTTIl proved, based on results in [3 1 and 1121 . that in the periodic case, 
the limiting particle distribution does not satisfy a linear Boltzmann equation. Caglioti 
and Golse GHEE) and Marklof and Strombergsson H3] QU \T7\ CDD independently found 
sharp estimates of the distribution of free path lengths and that the Boltzmann-Grad limit 
corresponds to a Boltzmann like equation in an extended phases space that does describe 
the evolution of a particle density in the limit of small s; whereas the results in j6l|2llH] 
are restricted to two space dimensions and rely on an independence assumption, fl5\ [161 
[l7l[T8ll provides a complete proof valid for any space dimension. Related results valid for 
in the two-dimensional case, can be found in J2][T), and in (4). Situations where scatterers 
are randomly place on a periodic lattice have been considered in |5 1 and |20|. 

Here we are interested in the case where the obstacles are distributed as the atoms in 
a quasi crystal. By definition a crystal is "a solid with an essentially discrete diffraction 
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pattern" ETTl . It has been known for a long time that periodic crystals in three dimen- 
sions must belong to one of fourteen symmetry classes. A quasi crystal is a crystal whose 
diffraction pattern exhibits a forbidden symmetry. The first reports on experimental results 
indicating that such solids exist were treated with suspicion by the scientific community, 
but in 201 1 their discovery was awarded the Nobel Prize in Chemistry [ 13 1. There are also 
many mathematical abstract constructions that give the same results, starting e.g. from the 
Penrose tiling (see e.g. ED). It is then a natural question to ask whether the Boltzmann- 
Grad limit of a Lorentz gas in a quasi crystal behaves more like the random or periodic 
case. In this paper we present simulation results on a one-dimensional model, which give 
strong support for the latter: the free path length distribution decays polynomially, just as in 
the periodic case, whereas in the random case, the path length distribution is exponentially 
decaying. 

The next section gives precise definitions of the Lorentz model, and of its Boltzmann- 
Grad limit, and some results concerning the periodic and random cases are discussed. 
Thereafter, in Section [3] a construction of quasi crystals is given; that section is mainly 
based on ETI . 

Section[4]presents the simulation method and the results, and the paper ends with con- 
clusions and some prospects for future research in Section[5] 

2. Free path length distributions in the Lorentz model 

The discussion in this section is restricted to the Lorentz gas in two dimensions, and 
hence we consider a point distribution r e C M 2 , or more precisely, a family of point 
distributions parametrized by e. At each point p E T s we put a circular obstacle with radius 
e and center at p, and then we study the motion of a point particle moving with constant 
speed, |v| = 1, along straight lines between the obstacles, and specularly reflected when 
hitting an obstacle. Hence the phase space for one point particle is f2p e x S 1 , where 

^r e = K 2 \ IJ B p (e) 
per 

and B p (e) is the closed ball of radius e centered at p g R 2 . For any initial point (xo, Vq) 
we define 

Zr.(*o,«o) = {x(t),v(t)), 
the position of the point particle at time t, taking into account all reflections on the set 
of obstacles. As long as the obstacles do not overlap, this is well defined for all t > 0. 
Moreover, the the free path length is defined as 

T e (x , v ) = inf {* > x + tv fip e } ■ 
We also define the path length distribution 

(1) M ]aM) = Hm W{(x,t>) e (nr.n go(JZ))x gMr 6 (x,,)E]a,6D} 
J reVJ ' U fl^oo m(n rc nB (R)) 

Here m(-) is the Lebesgue measure restricted to f^r e xS 1 . In cases where, at least formally, 
one can prove that for all intervals [a, b], <j> E (]a, b[) remains bounded from above and below 
when e — > one speaks of a Boltzmann-Grad limit. 

The two typical examples of point distributions T e considered in connection with the 
Lorentz gas are the periodic distributions T e = a(e)Z 2 , and a Poissonean random distribu- 
tion with intensity a(e). 

For a periodic distribution it is more natural to define the free path length distribution 
by restricting to a lattice unit cell rather than to a ball of radius R as in ([T]). Of course, in 
the limit R — > oo the result is the same. In this case the Boltzmann-Grad limit corresponds 
to choosing a(e) = e 1 / 2 . The path length distribution <p £ in the Boltzmann-Grad limit has 
been studied in J3][T2], and then, with very sharp bounds in (6]|2l[8|] using the theory of 
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Figure 1 . Computing the free path length in a periodic Lorentz gas is 
equivalent to to a discrete map. 

continued fractions, and flTl [161 [171 [18 1 using methods based on Ratner's theorem. The 
result is that <fr E ([T, oof) — >• C T _1 asymptotically for large T when e 0. 

The distribution T e is a Poisson distribution with intensity a(e) if and only if for any set 

A C M 2 , 

Prob ( no of points inA = k)= (m(A) " (£)) " e -™(A)a(e) 

and in also, for 4 C K 2 and B C M 2 with A {J B = 0, the number of points in T E [j A 
and r e (J B are independent random variables. Here the Boltzmann-Grad limit is achieved 
by setting a(e) — e^ 1 , and letting e go to zero. This distribution differs from the periodic 
one in several fundamental aspects. First, obstacles may overlap, and hence the trajecto- 
ries (x(t), v(t)) cannot always be continued uniquely. However, the measure of such, bad 
trajectories goes to zero when e 0. Secondly, before taking the limit R — > oo in the 
definition of the path length distribution, the expression in the right hand side of equa- 
tion (JTJ is a random variable, but one that converges to a deterministic value both when 
R — > oo, and when e — > 0. The path length distribution <fi e ([T, oof) converges to a dis- 
tribution <f>([T, oof) ~ exp(— cT). The Lorentz gas with Poisson distributed obstacles has 
been studied in detail by Gallavotti [9], who proved that if f £ ,o(x, v) £ L 1 (flr c x S 1 ) are 
densities of initial points for point particles, and f e (x, v,t) — E \f e fi(T^ (x, v))] , then, 
assuming that / e ->/o£ i 1 (M 2 x S 1 ), it follows that f E (x,v,t) — > f(x,v,t), which 
solves a linear Boltzmann equation. 

A different class of random distributions of scatterers can be constructed starting from 
a periodic distribution, by removing obstacles randomly, independently, with some prob- 
ability p(e). Although for a given, positive e, the behavior is rather different from the 
Poissonean case, this difference disappears in the limit as e — > 0, and in particular, it is 
possible to rigorously derive a linear Boltzmann equation starting from such distributions 
(see EE)]]). 

The following section gives examples of mathematical constructions of quasi crystals, 
and the corresponding Lorentz gas. However, computing long point particle trajectories 
in a quasi crystal Lorentz gas is very time consuming, and therefore the simulation results 
presented in this paper are performed on a one-dimensional discrete model. In the two- 
dimensional periodic case, the path length distribution can be computed almost exactly 
by analyzing the discrete map given by the consecutive intersections of a trajectory with 
lines parallel to the lattice containing lattice points, as in Figure [T] We assume here that 
a = 1. For an initial point (xo,v) with xq = (qo,q"o) sitting on a horizontal line as in 
the figure, and with v E S 2 having an angle 9 to the vertical line, the distance between 
two consecutive points is tan(0). To find the free path length of a trajectory in a periodic 
Lorentz gas is then (almost) equivalent to computing 

(2) k = min { j | dist(q j , Z) < y/s/2} . 

In fact, we will have r e (xa, v) = y/e(k ±1)/ cos (9. 
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3. Quasi crystals 



The material in this section is mostly taken from M. Senechal's book "Quasi Crystals 
and Geometry" [21 1, which gives a good introduction to the subject. 

There are several ways of construction point sets in K™ that satisfy the requirement for 
being a quasi crystal: "an essentially discrete diffraction pattern exhibiting a forbidden 
symmetry". One of these is the so called projection method: 

• Let CP be a point lattice in R k (usually the standard lattice), and let K be the unit 
cell in the lattice that contains the origin. 

• Let £ be an n -dimensional subspace of M. k , such that £ n C p — {0}, and let £ ± be 
the orthogonal complement of £. 

• Let II and H 1 - be the orthogonal projections on £ and £ ± respectively. 

• Let X — C p H (IT 1 - (If) © £ J, i.e. the set of lattice points contained in the strip 
(or cylinder) obtained by translating the unit cell K along £. 

The quasi crystal is A = Ti(X), the orthogonal projections of X on £. This set is dis- 
crete (inf { \pi — p%\ \ p\-,P2 G A} = ro > 0), non periodic and relatively dense (meaning 
that there is Ro > such that every sphere in £ of radius greater than Ro contains at least 
one point of A) (see Figure [2]). That A is discrete and relatively dense implies in particular 
that it is possible to define a scaling that corresponds to the Boltzmann-Grad limit. 

The theory of quasi crystals started with the (experimental) discovery of a material ex- 
hibiting a five-fold symmetry. Mathematically this can be constructed using the projection 



method, starting from the regular lattice in 
Q5 with 32 vertices 
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A unit cell in this lattice is the hyper cube 



6 {-1,1} 



where e% = (1, 0, 0, 0, 0), es = (0, 0, 0, 0, 1). A five fold symmetry is given by cyclic 
permutation of the unit vectors ej, a transformation that can be represented by the rotation 
matrix 
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which by an orthogonal change of variables becomes 
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In these new coordinates, £ will be the subspace spanned by (1, 0, 0, 0, 0) and (0, 1, 0, 0, 0). 
The quasi crystal will then be the the set 

n(£ p u (n- L (g 5 )©£)) c£. 

A one-dimensional quasi crystal can be defined in the similarly as a projection on a 
one dimensional subspace of M 2 . Particular examples, which are convenient for numerical 
computations, are the point sequences known as Fibonacci sequences. They are defined by 
taking 

£ = {x G M 2 \x ■ lo = 0} , 

where 

r ■, \ -u 1 + ^5 
lo = ( — 1, t) with r = . 
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Figure 2. A two dimensional representation of the projection method 
for constructing quasi crystals 



Let v = \lu I = Vl + r 2 . There is an explicit formula for computing the sequence of points 
that constitute the one-dimensional quasi crystal T obtained by this construction: 

m 

TV 



(3) T = II {x m } with x m = — H 

II Til 



where || • || denotes the distance to the nearest integer. This sequence has many interesting 
properties, and in particular it satisfies the criteria that defines a quasi crystal. A particular- 
ity is that the interval between two consecutive points are of exactly two kinds, short and 
long: 

1 I 1 

•Em •Em—l — 01 -r . 

v v VT 

As we have seen above, there is a very direct connection between the periodic Lorentz 
gas in two dimensions and a discrete model, at least when it comes to computing the 
free path-length distribution, illustrated in Figure [T] and Equation (|2]). Here we define the 
rescaled free path length 

(4) T e (q , v)=e min { j | dist(q + jv, F) < e/2} . 

(Note that compared to Equation Q the scaling factor y^i is replaced by e, the only reason 
being to make notation more clean). 

In the following section we simulate trajectories in order to estimate the path length 
distribution and compare with simulation results when T is replaced by oT lr L, and with a 
Poisson stream of points Xj with intensity a. The factor a is chosen so that all three point 
distributions have the same density, asymptotically: 

a = lim = — . 

-R-i-oo 2R v 

This is not equivalent to the free path length distribution in a two-dimensional Lorentz gas 
as it is in the periodic case, however. 

4. Simulation method and results 

Here we study the distribution of free path lengths of a point particle jumping on the real 
line, i.e. the number of jumps needed before the particle falls into an interval of width e 
with center at either the quasi crystal T, as defined previously, the periodic lattice ^Z, or 
on a Poisson distributed set of points, Y e with intensity r 2 jv. The intensity of the Poisson 
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distribution and the scaling of the periodic lattice are chosen in order that all three obstacle 
distributions have the same density. 

The random starting points q are chosen randomly, uniformly in the interval 
10000^/r 2 . Obviously, for computing the path length distribution in the random obstacle 
distribution this is not necessary, qo = would give exactly the same result. Moreover, in 
this case the distribution can be computed analytically in the limit of small e. That case is 
included in the simulation only for comparison. In the periodic distribution of scatterers it 
would be more natural to chose a random initial point uniformly in the interval [0, v/t 2 \, 
but by definition the quasi crystal is not periodic, and therefore it is natural to chose a larger 
interval. 

For the simulations presented here, v is chosen uniformly in [0, v/t 2 ]. To obtain a result 
completely equivalent to the periodic Lorentz gas in two dimensions, one should have 
taken v — -ps tan 9 with 8 uniformly chosen in the interval [0, 7r/4], but simulations with 
different jump length distributions show give very similar results, and this is not presented 
here. However, we have only tried distributions with bounded densities. 

The position of the point particle after j jumps is denoted qj = qo + vj, and the points 
x m G J- in the Fibonacci sequence are computed with the formula ([3]). The points in the 
Poisson distribution are computed independently for each trajectory q m . 

The simulation procedure is then 

(1) Chose qo and v randomly 

(2) Compute kjr e — min { j | dist(go + jv, J-) < e/2} 

(3) Compute fcz, e = min {j | dist(g + jv, ^yZ) < e/2} 

(4) Compute fcr\ e = mm {j | dist(qo + jv, T) < e/2}, where T is the random obsta- 
cle distribution. Note that for a given trajectory the free path time is a random 
variable, depending on the realization of the obstacle distribution. Only one real- 
ization of the random obstacle has been chosen. 

(5) Make a record of each of these path lengths, and repeat a large number N times. 

For all simulation results presented below, N = 2x 10 7 , and the results results have been 
used to estimate Prob[fc*. e > K], where * represents the different obstacle distributions. In 
connection with the Lorentz equation, the relevant quantity is the path length distribution 
expressed in time units, Prob[efc* i£ > T], which is estimated as 

Number of trajectories with fe* E > T e^ 1 

(5) (M[T,oo[) = ' 

The first example, shown in Figure [3] compares the free path length distributions with 
the three different distribution of scatterers that we consider here for e = 10~ 5 . The graph, 
plotted in logarithmic scale, show the 1/T-behavior of 0*, £ ([T, oo[) both for the quasi 
crystal T and for the periodic distribution, while the random distribution gives a different 
result (an exponential decay, as expected). 

Figure|4]shows the same as Figure[3]but with e = 1CP 4 (left) and e = 10~ 3 (right). Note 
that the path length distribution for the random obstacle distribution is exponentially only 
in of the exponential distribution for large T. This corresponds the fraction of trajectories 
with a jump length v < e which are stopped by the first obstacle they reach. 

Finally, Figure [5] shows the path length distribution cf>jr E ( [T, oo [) for the quasi crystal 
with three different values of e. The graph on the left shows distribution of the number 
of free steps , kjr e , and the graph on the right shows the three curves expressed as a 
function of T, as in Equation [5] That these curves almost coincide is a strong indication 
that cf>jr }e [T, oo [) converges to some distribution <j)j^[T, oo[) as e — > oo, just as in the 
periodic case. 
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Figure 3. The free path length distribution based on 10 7 trajectories 
and e = 1CP 5 . The curves show the result for the Fibonacci distribu- 
tion of scatters (green), the periodic distribution (blue) and the random 
distribution (magenta). 




FIGURE 4. The free path length distribution computed as in Figure[3]but 
with £ = 1CT 4 (left) and e = 1CT 3 (right) 

5. Conclusions 

The simulation results show that at least the one-dimensional discrete time Lorentz gas 
in a quasi crystal behaves essentially as the periodic Lorentz gas, and at least that the path 
length distribution is not exponentially decaying, as in the random case. The simulation 
examples all refer to the specific case of a Fibonacci sequence based quasicrystal. There 
are other constructions of onedimensional quasicrystals, and in particular one can construct 
a whole family of point distributions constructed as non-peridic sequences of two different 
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FIGURE 5. The free path length distribution in the quasi crystal, for 
e = 10~ 5 , 10~ 4 and 10~ 3 . To the left the curves represent the number 
of steps and on the right, the curves represent the scaled time before a 
trajectory hits an obstacle. 

intervals (the book by Senechal |21 1 presents some, and give many references). In partic- 
ular I wanted to see whether the special choice of the golden ratio as ration between long 
and short intervals would give a qualitatively different result from other, say trancendent, 
ratios. But the qualitative behaviour seems to be the same, and hence no simulation results 
are presented here. 

However, the simulations have been carried out in a much simplified model of the gas, 
and therefore it is not obvious that a simulation of a real (two-dimensional) Lorentz gas in 
a quasi crystal would give the same result. 

Because the quasi crystals considered here are constructed as projections of a regular 
lattice in higher dimension, it is possible that the methods of Marklof and Strombergsson 
referred to above would also give results in this case. This possibility is mentioned in lfT4ll . 
but as of now, no further results have been published. 

On the other hand, there are rather explicit results for a discrete Schrtidinger equation 
in a one-dimensional quasi crystal like the one studied numerically here 1 19 1. 
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